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The Kosambi-Cartan-Chern (KCC) theory represents a powerful mathematical method for the 
investigation of the properties of dynamical systems. The KCC theory introduces a geometric 
description of the time evolution of a dynamical system, with the solution curves of the dynamical 
system described by methods inspired by the theory of geodesics in a Finsler spaces. The evolution 
of a dynamical system is geometrized by introducing a non-linear connection, which allows the 
construction of the KCC covariant derivative, and of the deviation curvature tensor. In the KCC 
theory the properties of any dynamical system are described in terms of five geometrical invariants, 
with the second one giving the Jacobi stability of the system. Usually, the KCC theory is formulated 
by reducing the dynamical evolution equations to a set of second order differential equations. In 
the present paper we introduce and develop the KCC approach for dynamical systems described by 
systems of arbitrary n-dimensional first order differential equations. We investigate in detail the 
properties of the n-dimensional autonomous dynamical systems, as well as the relationship between 
the linear stability and the Jacobi stability. As a main result we find that only even-dimensional 
dynamical systems can exhibit both Jacobi stability and instability behaviors, while odd-dimensional 
dynamical systems are always Jacobi unstable, no matter their Lyapunov stability. As applications 
of the developed formalism we consider the geometrization and the study of the Jacobi stability of 
the complex dynamical networks, and of the A-Cold Dark Matter cosmological models, respectively. 
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I. INTRODUCTION 

Continuously time evolving dynamical systems are one of the fundamental theoretical tools for the mathematical 
modelling of the evolution of natural phenomena. They are widely used in every branch of natural sciences, and 
they provide extremely useful and powerful descriptions of the evolution of various natural processes. We dehne 
a deterministic continuous dynamical systems as a set of formal mathematical rules that describe the evolution of 
points (elements) included in a set S with respect to an external continuous time parameter t S T [l|. Even that 
generally one can take the time parameter as being discrete, or a continuous function [l|, in the following we restrict 
our analysis to the case of a continuous time parameter. More exactly, from a formal mathematical point of view, we 
define a solution of a dynamical system as a map [I| 

4> : T X S S, {t, x) 1 -^ 4>{t, x), (1) 

satisfying the condition •) o (/)(s, •) = cj)(t s, •), Vt, s G T. Usually, in the case of realistic dynamical systems that 
may be required to model natural phenomena, one must add some additional mathematical conditions to the above 
definition. 

In most of the applications in natural sciences the usefulness of the dynamical systems is largely determined by their 
predictive capabilities. On the other hand, the predictive potential of the dynamical systems is determined essentially 
by the stability of their solutions. We must also note that in realistic scientific investigations some uncertainties 
in the measured initial conditions in physical, chemical or biological systems do always appear, and they must be 
taken into account when formulating the basic mathematical equations of the considered processes. Therefore, in 
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order for a dynamical system to be useful in practical applications, it must provide essential information, and control, 
not only on the time evolution the system, but also on the possible deviations of the trajectories of the dynamical 
system with respect to a given trajectory, taken as reference. It is important to point out that from the study of 
a dynamical system we must obtain two essential information: (i) the local description of the stability, and (ii) the 
global description of the late-time deviations of the trajectories. 

The global stability of the solutions of the dynamical systems is described, by using a rigorous mathematical 
approach, by the well established and intensely investigated theory of the Lyapunov stability. In this standard 
approach to the problem of stability, the fundamental mathematical quantities of interest are the Lyapunov exponents. 
The Lyapunov exponents measure the exponential deviations of the trajectories of a dynamical system from a given 
reference trajectory ii. However, in trying to determine the Lyapunov exponents one should be aware that usually 
it is very difficult to obtain them in an exact analytic form. Therefore, in order to estimate the Lyapunov exponents 
one must use numerical methods. Presently many numerical methods have been developed for their computation, 
and they are used extensively in the analysis of the time behavior of the dynamical systems i, i . 

The methods of the Lyapunov stability analysis are well established, and they offer a clear understanding of some 
stability properties of the dynamical systems. However, in order to give a deeper insight into behavior of natural 
systems, alternative methods for the study of the stability must also be considered, and developed. Then, once such a 
new stability method is obtained, one can compare in detail the predictions and results of the new proposed method 
with the results obtained by using the Lyapunov linear stability analysis of the given dynamical system. 

One of the approaches that could offer important new insights for the study of the stability of the dynamical systems 
is represented by what one may call the geometro-dynamical analysis. One of the first examples of such an alternative 
study of the stability is the Kosambi-Cartan-Chern (KCC) theory. The ea^ formulation of the KCC theory was 
introduced in the important works of Kosambi Cartan Q and Chern Q, respectively. From a mathematical 
point of view the KCC theory is inspired by the geometry of the Finsler spaces, which represents its theoretical 
foundation. The KCC theory is based on the fundamental assumption that there is a one to one correspondence 
between a second order, autonomous or non-autonomous, dynamical system, and the geodesic equations in a Finsler 
space that can always be associated to the given system (for a recent review of the KCC theory see i). Essentially, 
the KCC theory represents a differential geometric approach to the variational differential equations describing the 
deviations of the entire trajectory of a dynamical system with respect to the nearby ones [^. In the geometrical 
description introduced by the KCC theory one can associate to each dynamical system two connections, the first one 
being a non-linear connection, while the second one is a Berwald type connection. Five geometrical invariants can be 
constructed with the help of the non-linear and Berwald connections, respectively. The most important of these five 
geometrical invariants is the second invariant, called the curvature deviation tensor. Its importance from the point 
of view of scientific and mathematical applications relies on the fact that it gives the so-called Jacobi stability of the 
given dynamical system [9Hl2|. Different physical, biochemical, or engineering systems have been extensively studied 
with the use of the KCCtoy PQ. 

The study of the two-dimensional autonomous systems, and of their stability properties, represents an important 
field of application of the KCC theory. In P, P| two dimensional systems of the form 

du,.dv, , , 

— = f{u,v), —=g[u,v), (2) 

were studied under the assumption that the point (0,0) is a fixed point, i.e. /(0,0) = g(0,0) = 0. In the approach 
introduced in PP, after relabeling v as x, and g{u,v) as y, and by assuming that gii|(o,o) 0, it follows that one 
can eliminate the variable u. Since (u, v) = (0, 0) is a fixed point, it follows from the Theorem of Implicit Functions 
that in the vicinity of (a;, y) = (0,0) the equation g{u, x) — y = 0 has a solution u = u{x, y) . Since x = g = gu f + gvV, 
an autonomous one-dimensional second order equation is obtained, which is equivalent to the system and which 
is given by 

x^ + g^{x,y) =0, (3) 

where 

g^{x,y) = -gu{u{x,y),x) f{u{x,y),x) - gy{u{x,y),x)y. (4) 

Hence the Jacobi stability properties of Eq. ([3|) can be studied in detail by using the KCC theory PP. Con¬ 

sequently, the KCC stability properties of the first order system dJ) can also be obtained, and the in the in depth 
comparison between the Jacobi and Lyapunov stability properties of the two dimensional dynamical systems can be 
easily performed. 

Usually, the KCC theory is applied to dynamical systems by reducing the evolution equations to a set of second 
order differential equations. An alternative view of the KCC theory, in which the theory is applied to a first order 
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dynamical system, was introduced in [24| . The properties of the two dimensional autonomous dynamical systems were 
investigated in detail, and the relationship between the linear stability and the Jacobi stability was also established. 
As a physical application of the proposed approach the KCC type geometrization of Hamiltonian systems with one 
degree of freedom was considered, and the KCC stability properties of Hamiltonian systems were discussed. 

It is the goal of the present paper to generalize and extend the perspective introduced in [2^ for arbitrary n- 
dimensional first order dynamical systems. By starting with an n-dimensional autonomous first order dynamical 
system, we transform it ’’directly” to an equivalent second order autonomous system of ordinary differential equations. 
This transformation is done by simply taking the time derivative of all equations of the dynamical system with respect 
to the time parameter. The obtained second order equations can be investigated with the use of the KCC theory, 
which allows a full geometric description of the system. In particular, the expressions of the non-linear connection, of 
the Berwald connection, and of the deviation curvature tensor are explicitly obtained. Hence an in depth comparison 
between the Lyapunov stability and the KCC stability of the critical points can be performed. In particular, we 
obtain the important result that only even dimensional first order dynamical systems exhibit Jacobi stability. We 
also discuss in detail the case of the three-dimensional dynamical systems, which are always Jacobi unstable, no 
matter their Lyapunov stability. 

As a hrst application of the obtained results we briefly point out the possibility of investigating the KCC stability of 
complex dynamic networks. The theory and applications of complex networks represents a very active field of research 
in science and technology. In particular, collective motions of coupled dynamical networks, such as stabilization, 
synchronization, robustness, diffusion etc., have received a great deal of attention recently 0, [26l - l^ . We show that 
a KCC type geometric description of complex network dynamics is always possible, and that this description allows 
the study of the Jacobi stability of the network. 

The theoretical results obtained in the study of the dynamical systems have also found important applications in 
the held of cosmology [29l - l3^ . The Lyapunov stability properties of some cosmological models have been considered 
in detail [3^ . The investigation of the critical points in the phase space for given cosmological models can provide 
very important information on the overall cosmological dynamics. The existence of cosmological attractors can be 
related to cosmological solutions indicating the early or late-time behavior of the Universe. The saddle equilibrium 
points may give some indications on the transient cosmological solutions [s^ . Hence, the study of the properties of the 
equilibrium points of the cosmological models can offer essential information on the model behavior. As a further step 
in the study of the cosmological models from a dynamical systems perspective we consider the KCC geometrization 
and the Jacobi stability properties of the ACold Dark Matter (ACDM) models. The KCC approach provides the 
possibility of a full geometric description of the ACDM cosmologies, and the study of their Jacobi stability properties. 

The present paper is organized as follows. In Section |H] we briefly review the Lyapunov stability of systems of 
ordinary differential equations. A brief presentation of the basic concepts and results of the KCC theory is presented 
in Section Hm The comparison of the Lypunov and Jacobi stability for n-dimensional dynamical systems is performed 
in Section El The case of the three-dimensional dynamical systems is discussed in Section m The applications of 
the KCC theory for the geometrization and stability study of complex networks are presented in Section I VII The 
possibility of the KCC geometrization of the cosmological dynamics, and the study of the Jacobi stability of the 
ACDM model is considered in Section IVHI We discuss and conclude our results in Section IVHII In the Appendix we 
briefly discuss the relation between the stability of the critical points located at zero, and the arbitrary critical points. 


II. LYAPUNOV STABILITY OF SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS 

In the present Section we briefly review following some results on the Lyapunov stability of the dynamical 
systems, described by systems of Ordinary Differential Equations (ODEs). We first recall that, in general, the 
stability of a system of ODEs is determined by the roots of its characteristic polynomial. Indeed, let us consider the 
system of ODEs 


dx^ 

dt 

dx"^ 

dt 

dx^-^ 

dt 

dx^ 

dt 


f {x^,x'^,...,x^) , 
f {x^,x^,...,x'^) , 

.5 

/” {x^,x'^,...,x'^) , 


( 5 ) 
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where /i, / 2 ,..., are smooth functions. Let us linearize the above system about a steady state (a:J, ojq, ccq ) by 
associating the linear system 


/ dx^ \ 

dL 


dx‘^ 

dt 


/x^ 


= A 


( 6 ) 


where A = J /"■) ^2 is the Jacobian matrix of the system ([5]), evaluated at the steady state, 


A = 


a/1 

a/1 

dfl 

dx^ 

dx'-^ ' 

dx^ 

dr 

a/" 

or 

dx^ 

dx'^ ' 

dx^ 


/' yytl /Tt2 'T*^! 


The solutions of Eq. m are given by 


/x^\ 

x^ 

/Cie^i‘\ 

C 2 e^ 2 t 

\x"/ 



( 7 ) 


where is a constant vector and (Ai, A 2 ,A„) are the eigenvalues of the matrix A, obtained as the 

roots of the characteristic polynomial p{X) given by 


det {A — XIn) = 0, 


( 8 ) 


where In is the identity matrix . 

Definition. A solution of the dynamical system (0) is called stable if and only if all the roots Xi, X 2 , ■■■, A„, of the 
characteristic polynomial p{X) lie in the left hand side of the complex plane, i. e., Re A < 0 for all roots X. 

If this is true, then x^{t) = —>■ 0 exponentially as t —>■ 00 for all i, and hence [x^,x^, ...,a:") = (0,0, ...,0) is 

stable to small (linear) perturbations. 

In the n-dimensional case, the characteristic polynomial reads 

p(A) = A” + aiA”-i + ... + a„ = 0, (9) 

where the coefficients at, i = I,2,...,n are real numbers. Without any loose of generality we can always assume 
a„ ^ 0, because otherwise we have A = 0 and we obtain a characteristic polynomial of order n — 1 with the zero order 
coefficient non-vanishing. 

The necessary and sufficient conditions for p (A) to have all solutions such that ReA < 0 can be written as 



Oi 03 
1 02 


Oi 

03 

05 

o„ > 0, Di — oi > 0, D2 — 

> 0, £>3 = 

1 

0 

02 

04 



Ol 

03 


( 10 ) 


Dk 


Ol 

03 


1 

02 

04 

0 

Ol 

03 

0 

1 

02 


0 0 ... Ok 


> 0 , 


( 11 ) 


for all k = 1, 2,..., n. 

Another method to study the stability problem is to consider the relations between the non-vanishing solutions of 
p (A). We have 


n 

S .— ^ ^ — Q-i, 

2=1 


(12) 
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XiXj — 02, (13) 

i,3,i=Xi 

(14) 

p := Ai...A„ = (-!)"■ a„. (15) 

By taking into account the values of these coefficients we can obtain important information about the stability of 
the system JS]). 

Remark (Descartes’ Rule of Signs) [s^ 

Consider the polynomial ([9]) with On > 0. Let m be the number of sign changes in the sequence of coefficients 
...,ao}, ignoring any coefficients which are zero. Then there are at most m roots of p(A) which are real 
and positive, and further there m, m — 2, m — A, ..., real positive roots [s^. By setting uj := —A this property gives 
important information about the possible real negative roots, which is essential for the study of stability. 


III. KCC THEORY AND JACOBI STABILITY OF DYNAMICAL SYSTEMS 

In the present Section we briefly outline the fundamental concepts, basic assumptions, and results of the KCC 
theory. We also introduce the relevant notations, and present the definitions of the relevant geometric objects (for a 
detailed presentation see Q and [i3)- 


A. The geometry of arbitrary dynamical systems 


In the following we assume that the dynamical variable a;*, i = l,2,...,n are defined on a real, smooth n-dimensional 
manifold M.. The tangent bundle of M. is denoted by TAl. Usually A1 is thought as i?", M. = i?", and consequently 
TM = TR^ = ii”. Let’s consider a subset U of the Euclidian {2n + 1) dimensional space i?" x i?" x . On 
ft we introduce a 2n + 1 dimensional coordinate system (a;b?/*,t), i = l,2,...,n, where (a;*) = (a;^,a;^,...,a:”), 
(y*) = and t is the usual time coordinate. We define the coordinates y* as 


/ dx^ dx'^ da;” A 

\dt'~dt'""~dr) ■ 


(16) 


A basic assumption in the KCC theory is that the time coordinate t is an absolute invariant. Therefore, within the 
framework of our approach the only admissible coordinate transformations are of the form 


t = f, i* = i* (a;^, a;^,..., a;") ,i G {l,2,...,u}. 


(17) 


In many situations of scientific interest the equations of motion of a dynamical system can be derived from a 
Lagrangian function L : TA4 ^ R via the Euler-Lagrange equations. 


d dL 
dt dy^ 


dx^ 


= Fi,i = 1, 2,..., n, 


(18) 


where in the case of mechanical systems E), j = 1, 2,..., n, represents the external force. If the Lagrangian L is regular, 
by means of simple computations one can show that the Euler-Lagrange equations Eq. m are equivalent to a system 
of second-order ordinary differential equations given by 


dt"^ 


+ 2G® {x^ 


y^,t) = 0,1 G {l,2,...,u}, 


(19) 


where each function G® [x^,y^,t) is C°° in a neighborhood of some initial conditions {xo,yo,to) in fl- 

The basic idea of the KCC theory is that if an arbitrary system of second-order differential equations of the form 
(HU) is given, with no a priori Lagrangean function known, we can still study the behavior of its trajectories by using 
differential geometric methods. This study can be performed by using the close analogy with the trajectories of the 
standard Euler-Lagrange system. 



B. The non-linear and Berwald connections, the curvature deviation tensor, and the KCC invariants 

associated to a dynamical system 
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In order to study of the geometrical properties associated to the dynamical system defined by Eqs. dm), on the 


base manifold A4 we introduce first the nonlinear connection N, 


with coefficients Nj 


defined as [25 1 



( 20 ) 


From a geometric point of view the nonlinear connection A)* can be interpreted in terms of a dynamical covariant 
derivative as follows. Given two vector fields v, w defined over a manifold A4, we introduce the covariant derivative 
as Q 


V,% = 


+ Nj{x,y)w^ 


d 

dx^ 


( 21 ) 


For Nil{x,y) = Tji(x)y’', from Eq. (1211) we immediately recover the definition of the covariant derivative for the 
special case of a standard linear connection, as usually defined in Riemmannian geometry. 

On the open subset 0 C i?" x i?" x on which the dynamical system is defined, and for the non-singular coordinate 
transformations defined by Eqs. (fT71) . we introduce the KCC-covariant differential of a vector field CM through the 
definition [IMl 


DC 

dt 


dt 


-mc- 


( 22 ) 


For C = y* we obtain 


Dv" 

= Ny -2G^ = -e\ (23) 

dt ■' 

e®, representing a contravariant vector field, defined on the subset 12 of the Euclidian space, is called the first KCC 
invariant. Erom a physical point of view, and in the framework of classical mechanics, the first KCC invariant e® may 
be interpreted as an external force. 

As a next step in our study we vary the trajectories a;® (t) of the system (1191) into nearby ones according to the rules 


X® (t) =x\t)+ 7 ]C{t),f (t) 


y"{t) + v 


dCjt) 

dt 


(24) 


where \ri\ « 1 is a small parameter, and C{i) ^re the components of a contravariant vector field, which is defined 
along the trajectory x®(t) of the dynamical system. By substituting Eqs. (1^ into Eqs. (ITOl) . and by taking the limit 
r? —>■ 0, we obtain the so-called deviation, or Jacobi, equations as given by (l0l - [l^ 


d'^C 

dC 


'"if+ = 


(25) 


By using the KCC-covariant derivative, defined by Eq. (1^ . we can rewrite Eq. (ESI) in a covariant form as 


D^C 

dC 


= py, 


where we have denoted 


aC® , ■ , dN] . , dm 


w y 

In Eq. (1221) we have also introduced the tensor G®;, defined as 

dm. 

/^l — _ 

- Qyl ’ 


(26) 


(27) 


(28) 


and which is called the Berwald connection. 
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The tensor is the second KCC-invariant, and it may be called alternatively as the deviation curvature tensor. 
We shall call Eq. (1^51) the Jacobi equation. In either Riemann or Finsler geometry, when the system of equations (fTOl) 
describes the geodesic motion, then Eq. (12611 is the Jacobi field equation corresponding to the given geometry. 

A scalar quantity which can be constructed from Pj is the trace P of the curvature deviation tensor, which can be 
obtained from the relation 




■ nin^ 


dNl 


(29) 


In the KCC theory one can also introduce the third, fourth and fifth invariants of the second order system of 
equations (HH). These invariants are defined according to 0 



dyi 


P* = 
? ^jkl — 


dP: 


jk 


dy’- 


n® = 

^jki - 


dG\ 


jk 


dy^ 


(30) 


From a geometrical point of view the third KCC invariant can be interpreted as a torsion tensor. The fourth 
and fifth KCC invariants and represent the Riemann-Christoffel curvature tensor, and the Douglas tensor, 
respectively I,111. It is important to point out that in a Berwald space these tensors always exist. In the KCC theory 
these five invariants are the basic mathematical quantities describing the geometrical properties and interpretation of 
an arbitrary system of second-order differential equations. 


C. Jacobi stability of dynamical systems 

In many scientific applications, including the study of the stability of physical, chemical, biological or engineering 
systems, the study of the behavior of the trajectories of the dynamical system given by Eqs. (IT^ in the vicinity of 
a point cc* (to) is of extreme importance. In the following for simplicity we take the origin of the time coordinate at 
to = 0. We interpret the trajectories = x^{t) of the system of equations (fTOl) as curves in the Euclidean space 
(P", (.,.)), where (.,.) represents the canonical inner product of P". As for the deviation vector ^ we assume that it 
obeys the set of initial conditions ^ (0) = O and ^{0) = W ^ O, respectively, where O G P" is the null vector [Ml- 

To describe the focusing/dispersing tendency of the trajectories around to = 0 we introduce the following mathe¬ 
matical description: if the condition ||^ (t)|| < t Ri O'*' is satisfied, then the trajectories are bunching together. On 
the other hand, if the deviation vector satisfies the condition ||^ (t)|| > t Ki 0+, it follows that the trajectories of 
the system m have a dispersing behavior [Ml. The focusing/dispersing behavior of the trajectories of a dynamical 
system can be also described from a geometric point of view by using the properties of the deviation curvature tensor 
as follows. For t fv 0+ the trajectories of the system of equations Eqs. (ITOl) are bunching/focusing together if and 
only if the real parts of the eigenvalues (characteristic values) of the deviation tensor P^ (0) are strictly negative. On 
the other hand, if and only if the real part of the eigenvalues (characteristic values) of P^ (0) are strictly positive, the 
trajectories are dispersing [Mi- 

Based on the above qualitative discussion, we introduce the rigorous definition of the concept of Jacobi stability for 
a general continuously time evolving dynamical system as follows [Ml: 

Definition: Let’s assume that with respect to the norm ||.|| induced by a positive definite inner product, the system 
of differential equations Eqs. m satisfies the initial conditions 

I la;* (to) - i* (to) 11 = 0,1 lx* (to) - x* (to) 11 ^ 0. 


Then the trajectories of the dynamical system given by Eqs. m are called Jacobi stable if and only if the real 
parts of the characteristic values of the deviation tensor P^ are strietly negative everywhere. If the real parts of 
the characteristic values of the deviation tensor Pf are strictly positive everywhere, the trajectories of the dynamical 
system are called Jacobi unstable. 


IV. RELATION OF JACOBI AND LYAPUNOV STABILITY FOR HIGHER DIMENSIONAL 

DYNAMICAL SYSTEMS 

In the following we will always assume, without loosing the generality, that the fixed point of the system of equations 
Eqs. ([5]) is (0,0,...,0) (see the Appendix for a discussion of the case of arbitrary fixed points). Let us consider the 
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system of ODEs given by Eqs. We point first out that, geometrically, we can describe the solutions of Eqs. 
as a flow ipt : D C i?" —!> i?", or, more generally, (ft '■ D C A4 ^ Ai, where is a smooth n-dimensional manifold. 
The canonical lift of (ft to the tangent space TAi is given by 

(ft ■ TAi TAi, (ft (u) = {(ft (u) , (ft (u)) , (31) 


where 


{u) 


df {t, u) 
dt 


(32) 


A. Geometrical description of higher dimensional dynamical systems 


In terms of dynamical systems we simply take the derivative of ([5]) with respect to the time parameter t, thus 
obtaining 

= /i (a:\...,a;’") + /j (a;\...,x") + ... + (x\...,x") y", 




dt^ 


= fi (x\...,x") + /I (x\...,x") y^ + ... + (x\...,x") y", 


where 


dt^ 


= A" {x\ ..., x") yi + A" {x\ ..., x") y^ + ... + (x^, ..., x") y", 


. df{x\...,x^) . . tdx^ 

Jj = -= 1 , 


In other words, on TAi we obtain 


dt 




2 = 1 


(33) 


dy"^ 

dt 


= ^/^(x^...,x")y^ 


2=1 


where (x^, x^,..., y^,..., y") are local coordinates on TAi. Hence one can see that the above system is actually a linear 
dynamical system on the fiber 3,2 Al. The system can be written as 


d^x* 

dt^ 


y'' = 0 ,i = 1 , ...,n. 




By comparing Eq. (IMl) to Eq. (ITOl) we find 

1 ” 

G* = --^A(xi,..,x”) y^^ = !,...,! 

or, equivalently, 


A;=l 


>02 


W) 






fy"\ 


\y^J 


= —-J • y. 
2 ^ 


(34) 


(35) 


(36) 


VELi/fe (a;i,..,x")yV 
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Using Eq. (EOl) we obtain 


or 


and hence 


m = - 
^ 2 


ldf{x\..,x^)_ 1 


dx^ 


,i,j = 


(^) 


3 / = 


'Ni Ni... Ni: 


,iVf ... Nil 




(37) 


(38) 


dm 

Then, for the components of the deviation curvature tensor (U*), given by Eq. (l27l) . we find 

P)r<i r)N'^. 

Then, with the use of Eqs. (I55t and dSZD we obtain 

. n 1 ^ 


(39) 


(40) 




1^1 


Therefore we have obtained the following 

Proposition 1. The curvature deviation tensor associated to an n-dimensional dynamical system is given by 


1 1 


(41) 


k,l^l 


^ fk fk £k 

Jll Jl2 ••• Jin 

where Hfk = I . is the Hessian of [x^,x^, ...,x^^. 

£k £k £k 

\.J nl J n2 J nn/ 


B. Lyapunov stability versus Jacobi stability 


Evaluating P- at the stability state (0, ...,0) we obtain 


(^)l 


( 0 ..... 0 ) 



+ E/E) 

( 0 ..... 0 ) 


( 0 ..... 0 ) 




(42) 


where A is the Jacobian matrix, with coefficients evaluated at the fixed point. We briefly discus the case when the 
stability state is not at zero in the Appendix. 

Let us now recall some general results from linear algebra. 

Lemma 1: Let A be an (n x n) matrix, and denote by Ai, A 2 ,..., A„ its eigenvalues. Then 

(i) The eigenvalues of kA are kXi,...,kXn, for any scalar k ^ 0. 

(ii) The eigenvalues of A^ := A x . . . x A are (Ai)^ ,..., (A„)^. 

k times 

Proof. Let us recall that if A is a complex matrix with complex eigenvalues Ai,..., A„, counted with their complex 
multiplicities, then always there exists an invertible matrix Q such that the matrix T := QAQ~^ is upper triangular, 
and the entries on the diagonal are exactly Ai,..., A„. This means that A and T have the same characteristic polynomial 
Pa = Pr = (A — Ai)... (A — A„) = 0 (see for instance the book [1^). Using this observation the proof of the Lemma 
is simple. Indeed, remark that T = QAQ~^ implies kT — Q (kA) Q~^, where fc is a scalar. Then, remarking that 
diagonal entries in kT are fcAi,..., fcA„, we get p^i = pr = (A — fcAi) (A — kX 2 )... (A — fcA„), and therefore (i) is proved. 
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Similarly, by remarking that T = QAQ ^ implies = Q ■ ■ Q and that the diagonal entries of are 

(Ai)^ ,(Afc)^, then (ii) follows by the same arguments as above. 

From Lemma 1 we obtain an important result that relates the two types of stability. 

Lemma 2 : If Ai, A 2 ,..., An are eigenvalues of the matrix A = J (/i, / 2 , fn) |(o,...,o) then fii = (Ai)^,..., = (A„)^ 

are eigenvalues of the deviation curvature matrix P|(o,...,o) = l(o,...,o)i hj = 1 ; 

This Lemma leads to a very natural relation between classical stability theory and KCC theory. 

Theorem 1. Let us consider the ODEs system 0), and denote by A := J (/i, / 2 ,..., fn) |(o,...,o) t/ie Jacobian matrix 
evaluated at the fixed point (0,..., 0). 

1. If A has only real eigenvalues, that is, the steady state is a (stable or unstable) node, then P|(o 0 ) has only real 
positive eigenvalues, that is, the steady state is always Jacobi unstable. 

2. If A has conjugate complex eigenvalues p,j := aj + ijlj, p, := aj — iflj for some j so that 2j < n, then 

2.1. If 2j = n, then the steady state Jacobi stable if and only if — (3'^ < 0, for all j = 1,..., n/2; 

2.2. If 2j < n, thenp(X) has real and complex eigenvalues and therefore the steady state is Jacobi unstable. 

We can now summarize our main results as follows. 

Theorem 2. The steady state (0, ...,0) of the dynamical system (0) is Jacobi stable if and only if satisfies the 
following 

(i) It is even dimensional, i.e. n = 2j. 

(ii) All eigenvalues of the characteristic polynomial p{X) = 0 are complex conjugate, i.e. Xj := aj + i(3j, Xj := 

j = 

(Hi) For any j = l,...,n/2 we have 


of - /3f < 0. 


V. THE CASE OE THE THREE VARIABLE DEPENDENT DYNAMICAL SYSTEM 


In the following we consider the dynamical system 


dx^ 

dt 


/(a 




dx"^ 

dt 




(43) 


— = h (x^x^x^) , 

and consider again that the steady state (xJ, Xq, Xq) = (0,0,0). The characteristic polynomial associated to the above 
system reads 

p (A) = X^ — s ■ X"^ + fj. ■ X — p, (44) 

where we have defined 

s := Ai + A 2 + A 3 , (45) 

p, := A 1 A 2 + A 2 A 3 + A 3 A 1 , (46) 


p — A 1 A 2 A 3 . 


(47) 


A simple computation of the determinant 


IA-AJ3I = 0, 


(48) 
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where 


j ail 

012 

Ol3 

A= 021 

O 22 

O 23 

Vosi 

O 32 

033 


^ df df[x^,x^^x^^ df(^x^,x^^x^^'^ 

Qj^l P»^2 


dx'^ 


dx'-^ 


dg(^x^ ,x^ ,x^^ dg(^x^ ,x^ ^x^^ dg(^x^ ,x^ ,x^'j 


“aP” 


dx'-^ 


dh(^x^ ,x^ ,x^^ dh(^x^ ,x^ ,x^^ dh(^x^ ,x^ ,x^^ 


dx'^ 


dx'^ 


dx^ 


( 0 , 0 , 0 ) 


and I 3 is the (3,3) identity matrix, implies the familiar form 


p (A) = — tr A • i (tr Af — tr (A^) 


• A — det A = 0, 


(49) 


(50) 


where tr A and det A are the trace and determinant of the matrix A, and tr is the trace of the matrix A^ := A-A. 
The classification of the steady states of a three-dimensional dynamical system was given by Poincare, and can be 
summarized as follows: 

1. Ai, A 2 , A 3 real roots of same sign of p (A): p is called a node. 

2. Ai, A 2 , A 3 real roots of p (A), but not of same sign: p is called a saddle. 

3. Ai, A 2 G C, A 3 S R, such that Re (Ai_ 2 ) and A 3 are of the same sign: p is called a focus. 

4. Ai, A 2 G C, A 3 G M, such that Re (Ai, 2 ) and A 3 are not of the same sign: p is called a saddle-focus. 

5. Ai, A 2 G C, A 3 G M, such that Re (Ai, 2 ) = 0: p is called a center. 

This classification can be combined with the Lyapunov stability in the case Re (A^) < 0 for i = 1,2,3. Taking now 
Theorem 2 into account it follows 

Theorem 3. The fixed point p of the three-dimensional dynamical system o is always classified as Jacobi 
unstable, regardless its Lyapunov stability. 

Remark. Our Theorem 3 above might be understood that in the case n = 3, and more generally n = 2j, j G N, 
KCC theory gives no useful information on the steady states. This is completely true. Without having any Jacobi 
stable trajectory, if p is a focus or a saddle-focus, then the case when Ai, A 2 G C, A 3 G R, < 0, where 

Ai ,2 = a ± i/3, gives information about a Jacobi-type saddle-focus, distinct from the usual Jacobi unstable behavior. 
We will not pursue here further investigations of these behaviors. 


VI. JACOBI STABILITY ANALYSIS OF DYNAMICAL NETWORKS 


The study of the dynamics of complex networks is has many applications in both science and engineering. Different 
processes like such as stabilization, synchronization, robustness, diffusion etc. can be described by using networks of 
coupled swtems. From mathematical point of view a complex network of coupled systems is described by the system 
of ODEs in 


dx^ 

dt 


F* (a 


(a 




) = 




),i = 1 , 2 , 


(51) 


where a:* with f G 1, 2,..., N are dynamical variables, F and iJ are evolution and coupling functions, respectively, a is 
a constant, and L* is the Laplacian matrix, defined by L- = F (the connectivity degree of node i), L\. = —1 if nodes 
i and j are connected, and L* = 0 otherwise. 

A Lyapunov linear stability analysis of the dynamical system describing networks of coupled systems can be per¬ 
formed by (i) expanding around a critical state Xi = X 2 = ■■■ = xn = with solution of = F (x^,.) = 0, (ii) 
diagonalizing L to find its N eigenvalues 0 = Ai < X 2 < ■■■ < ^n, and, finally, (iii) writing the governing equations 
for the normal modes yt of the perturbations 

y, = [F'{xl)-aAH'ixl)]y.. (52) 


The above equations have all the same form, but different effective couplings a = aXi. 

As a first step in the Jacobi stability analysis of the network dynamics we need to extend the dynamical process 
given by Eq. (IHTll to differentiable manifolds. By taking the derivative of Eqs. m with respect to the time we obtain 


d^x^ _ ^ I 
dt^ dx’- ^ 


dF^ 
dx^' 


-aLi 


dH^ 

dx^ 



dF^ 


(53) 












12 


Eq. (l53)) has the same mathematical form as Eq. dun, with 




1 /9F* 


dH^ 


dx'^ dx 


(54) 


Therefore we can apply now the KCC theory for the description of the dynamics of the network system. The non-linear 
connection associated to the system Eq. (I55|) is given by 


^ dy^ 


1 fdF^ 


-aLl 


dH^ 


'” = -2 


1 f dF^ 


dH^ 


o ■ - ctT* . 
dxi dxi 


(55) 


2 \ 9a;" " ’’ 

All the components of the Berwald connection G®; = dN^/dy’" are identically equal to zero. Then we easily obtain 


- 2 ^=r 


92 E® 


dx^ V dx^dx'^ 


-aLl 


d'^H^ 

dx^dx'^ 


dm 

y ^ 


1 / 9^f® _ 

2 V9a;'’9a;" ^ ''dx^dx^ 


giving for the deviation curvature tensor the expression 


1 / 92 ^® d'^m 

— aL^- 


P" = , 

^ 2 V dx^ 9a;" 


dx^dx^ 


1 /9E® , 9iJ 

9x* ^ ’’ 9x* y V^a;'^ 


fdF^ 


aLl 


dm 


The trace of the deviation curvature tensor of the complex network is given by 


P = P! = 


1 / 92 ^® 


2 V 9a;® 9a;" 


-aLl 


d^m 
9a;® 9a;" 


„ 1 /9F® ^,977’- 
^ ^ 4 ( 9x' 9x' 


9F' , 977® 

9a;® ^ ^ 9a;® 


(56) 


(57) 


(58) 


Therefore a dynamical network is Jacobi stable if and only if the real parts of the eigenvalues of the deviation 
tensor PJ are strictly negative everywhere, and Jacobi unstable otherwise. By assuming that the critical point of the 
network is given by (0,..., 0), it follows that the deviation curvature tensor evaluated at the critical point has the form 


ml _ 1 /9F® ,977®\/9F' ,977® 

^j l(o....,o) “ 4 dx^ ) dx3 


(59) 


(o,....o) 


Therefore all our previous results, summarized in Theorem 2, can be applied for the case of complex networks. 
In particular it follows that only even-dimensional networks are Jacobi stable. The geodesic deviation equation can 
also be obtained, and its study may give some insights into the development of chaos in complex networks. 


VII. KCC GEOMETRIZATION AND THE JACOBI STABILITY OE THE ACDM COSMOLOGICAL 

MODELS 


In the following we restrict our analysis to isotropic and homogeneous cosmological models, described by the 
Friedmann-Robertson-Walker line element, given by 


ds^ = —dt^ + a?{t) [dx^ + dy"^ -|- dz^^ , 


(60) 


where a{t) is the scale factor of the Universe. The cosmological expansion for a Universe filled with pressureless dust 
and radiation is described by the Friedmann equations, which take the form 


3772 — -|- pr + a, 77 


1 

2 




(61) 


Pm+iHpm =Q,pr +'^Hpr =Q (62) 

where a dot denotes the derivative with respect to the time t. In Eqs. (IHTll and (l62|) 77 = d/a is the Hubble function. 
Pm is the matter energy density, pr is the energy density of the radiation, and A is the cosmological constant. 
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A. Friedmann equations as an autonomous dynamical system 


In order to formulate the cosmological evolution in terms of a dynamical system, we introduce first the matter, 
radiation and cosmological constant density parameters (17,^, ^a)i defined as 


— 


Pn 


3i?2’ 




3iJ2 


, r^A = 


A 


The density parameters satisfy the relation 


+ fir + Ha — 1- 


(63) 


(64) 


As the basic variables {x,y) in the phase space we choose x = fir, and y = Ha, respectively. Then flm = \ — x — y, 
and the range of the variables is since 0<a:< l,0<j/< 1, and 0 < H^ < 1. The physically relevant phase space for 
the cosmological dynamics is defined as = {(x, y) : x + y < 1, 0 < a: < 1,0 < y < 1}. By taking the time derivatives 
of X and y with respect to the time, after introducing the new time variable r = lna(^, it follows that the Friedmann 
equations can be formulated in terms of an autonomous dynamical system given by 

dx 
dr 


-x(I - X -1- 3y), 

(65) 

= (3 -k X - 3y)y. 

(66) 


dr 

The critical points of the system and (IMl) in 'k can be obtained by solving the algebraic equations a;(l —a;+3y) = 0 
and (3 + a: — 3y)y = 0, respectively, and are given by 

PdS = {x = 0,y = l},Pr = {x = l,y = 0}, = {x = 0, y = 0}. (67) 

The Lyapunov stability properties of the system can be obtained from the study of the Jacobian matrix 


J = 


— 1 + 2x — 3y —3x 
y 3 + X — 6y 


( 68 ) 




The critical point P^js = (0,1), Ha = 1, with a(t) oc exp ^y'A/3tj, is associated to an accelerated de Sitter type 

expansion. It is a future attractor [s^. The critical point Pr = (1, 0), with H^ = 1 and a{t) oc v^, corresponds to the 
radiation-dominated era, and is a source point, or a past attractor. Finally, the critical point Pm(0,0), with Hm = 1 
and a{t) oc , corresponds to the matter-dominated phase of the cosmological expansion, and to a decelerating 
expansion. Pm is a saddle critical point [s^. 


B. The KCC geometrization and the Jacobi stability of the Friedmann equations 


In order to apply the KCC theory to the dynamical system given by equations (1651) and (1661) . we relabel the variables 
as X = x^ and y = x^. We also denote y^ = dx^ jdr and y^ = dx^/dr, respectively. Hence we obtain 


dx^ 

dr 


= —x^ (l — x^ -I- 3x^) = / (x^, x^) , 


(69) 


dx^ 

dr 


(3 -I- x^ — 3x^) x'^ = g (x^,x^). 


(70) 


Then, by taking the derivative with respect to r of Eqs. dMl) and dlO]), we obtain the following lift on the tangent 
bundle of the cosmological dynamical system. 


d^x^ 

dr^ 

d^x^ 

dr^ 


(—1 -1- 2x^ — 3x^) y^ — 3x^y^, 

(71) 

= x^y^ + (3 -1- x^ — 6x^) y^. 

(72) 







14 


By comparison with Eqs. (HU) we obtain immediately 

^ ) 2 x'^y^ + (3 + a;^ — 6x^) ) ’ 


(73) 


The components of the non-linear connection (TV) = Nf associated to the cosmological dynamical system are 
obtained as 


1 /—I -I- 2x^ — 3x‘^ —3a:^ 


{N)=N^ = -, 


For the components of the deviation tensor P = P/ we obtain 


3 -I- a;^ — 6 a;^ / 


P=^{Hfy H,-yY + ^J^f,g)., 


(74) 


(75) 


where Hf = 


fii fr. 
/21 / 2 ; 


is the Hessian of /, and Hg is the Hessian of g. By evaluating P at the critical points gives 


P\(y^=0,y^=0,x^=xl^,x^=x'^^) ’ 


(76) 


where A = J\(^yi^Q y 2 ^Q x'^=x^ )■ Explicitly, the curvature deviation tensor for the Friedmann cosmological 

dynamical system can be obtained as 


i(-2a;^ -I- 3a;^ -|- 1)^ -|- 


2 I 3x^x^ 


PI _| -TOO, -rxj -I- 2 — |a;^(a;^ — 6a;^-I-3)-I-|a:^(—2a:^-I-3a;^-I-1) 

l(yi=0,y==0) - I (^1 _ 6^2 + 3 ):r 2 + 1 (_2^1 + 3^2 ^ ^)^2 (^1 _ g ^2 ^ 3)2 ^ 3PP 


(77) 


We compute the numerical values of the curvature deviation tensor first at the critical point Pds = = 0, a;^ = l). 

In this critical point the curvature deviation tensor takes the form 


PI _ ( 4 0 

^ l(i/i= 0 .i/ 2 =o,a:i= 0 .a :2 = l) ~ 1 g 


(78) 


and has the eigenvalues Ai = 4 and A 2 = 9. Hence the critical point P^s of the standard ACDM cosmological model 
is Jacobi unstable. For the critical point Pr = [x^ = IjX^ = O), we obtain 


^\{y'^= 0 ,y'^= 0 ,x^ = l,x'^= 0 ) I Q 2*5 


1 21 
4 4 


(79) 


with the corresponding eigenvalues Ai = 16 and A 2 = 1/4. Hence the critical point Pr is also Jacobi unstable. Finally, 
for the critical point Pm = (0, 0) we obtain 


E*| (yl=0,1/^=0,a;i=0,a;2=0) y 0 9 J ’ 

with the eigenvalues Ai = 9, A 2 = 1/4. Hence all the critical points of the dynamical system associated to the 
Friedmann equations are Jacobi unstable. 


VIII. DISCUSSIONS AND FINAL REMARKS 


In the present paper we have extended and generalized an alternative approach to the standard KCC theory 


for first order autonomous dynamical systems, introduced in [2^. This approach is based on a different, ’’direct” 


transformation of the first order n-dimensional system to second order differential equations. We have presented and 
discussed this ’’direct” approach in detail for arbitrary dimensional dynamical systems. For this class of arbitrary 
dimensional dynamical systems we have developed and discussed in detail the two basic stability analysis methods - 
the (Lyapunov) linear stability analysis and the Jacobi stability analysis. From the point of view of the KCC theory, 
such a ’’direct” approach allows for a generalization of the geometric framework for first order systems. Consequently, 
the parameter space is increased, as well as the predictive power, of the method. While in the two dimensional case 
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there is a good correlation between the Lyapunov stability of the critical points, and the Jacobi stability of the same 
critical points [T^ . a very different picture emerges for dynamical systems with n> 2. The fundamental result of our 
analysis indicates that only even-dimensional dynamical systems are Jacobi stable, while odd-dimensional dynamical 
systems are always Jacobi unstable, no matter the Lyapunov stability of their critical points, describing the robustness 
of the corresponding trajectory to a small perturbation [l^ . On the other hand, it is important to point out that 
the Jacobi stability analysis is a very convenient way to describe the resistance of limit cycles to small arbitrary 
perturbation of trajectories. 

As a first application of our approach we have considered the possibility of a geometric description of the dynamics 
of complex networks. The dynamical behavior of such systems is described by n-dimensional dynamical systems, 
consisting of systems of first order, highly non-linear, differential equations. Such systems can be easily interpreted in 
the geometric framework of the KCC theory, and the associated geometric quantities (non-linear connection, Berwald 
connection, and deviation curvature tensor) can be obtained easily. The results obtained in this study allow a full 
geometric approach to the stability of complex systems, as well as the investigation of the Jacobi stability of the 
critical points. Hopefully, future studies may shed some light on the relation between Jacobi stability and chaotic 
behavior in complex networks. 

We have also considered the Jacobi stability analysis of the ACDM cosmological models. These models can be 
naturally represented in terms of an autonomous second degree dynamical system. The KCC theory allows a straight¬ 
forward geometrization of this system, which can be described in purely geometric terms with the help of a non-linear 
connection, the associated covariant derivative, and the curvature deviation tensor. The eigenvalues of the deviation 
curvature tensor as estimated at the critical points give the Jacobi stability properties of the cosmological model. 
We have considered three such models, corresponding to the de Sitter expansion, and to the matter and radiation 
dominated epochs, respectively. Our results show that the critical points corresponding to this cosmological models 
are all Jacobi unstable. 

In summary, in the present paper we have introduced and studied in detail some geometrical methods necessary for 
an in depth analysis and description of the KCC and Lyapunov stability properties of arbitrary dimensional dynamical 
systems. Such methods may help scientists to better understand the time evolution and the dynamical properties of 
natural phenomena. 
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Appendix A: The case of the arbitrary fixed points 

For the sake of simplicity we have assumed in our exposition that O = (0,..., 0) is the fixed point of ours system of 
ODEs. This assumption is not essential, and all the presented results remain true for arbitrary fixed points. Again 
for the sake of simplicity we will explain this fact only for the two-dimensional case. The higher dimensional case can 
be treated in an identical way. 

Let us consider the autonomous system of ODEs 


and consider the variable change 


du 

dt 




dv 

dt 


9 {u,v ), 


(Al) 


u := u — uo,v = V — Vo, (A2) 

where uq, vg are some real constants. For these new variables the system of ODEs Eqs. become 


— = f (u + uo,v + Vo) f {u,v), — = g{u + ug, V + Vg) =: g {u, v). (A3) 

It is easy to see that ODEs (IA3I) and are topologically and geometrically equivalent, hence they have the same 
stability behavior at the fixed points. 

Indeed, geometrically we recall that giving mil is equivalent in defining a vector field 


X\(^u,v) ■= fiu,v)—+ g{u,v) 


' dv' 


(A4) 


A solution of Eqs. (lAip is given by the flow (p^ of X in (or, more generally, in the two-dimensional manifold Ai. 
Similarly, (IA3I) gives a vector field X\(^u,v) obtained from X by parallel translation given by (IA2|) . Since du/dt = du/dt, 
dv/dt = dv/dt, it is clear that the direction of X does not change. 

The geometrical underlying reason for this result is clear. The group of isometries of the Euclidian plane is 
made of translations and rotations, and hence it is clear that the translation (IA2I) actually maps any solution curve 
of (lAll) to a solution curve of (IA3I) by keeping the orientation. 

From topological point of view, we recall that two flows ip* and on some manifold M are topologically equivalent 
if there exists a homomorphism h : j\4 —^ Ai that takes the trajectories of to the trajectories of if}*, while 
preserving the orientation. In our case, due to the linearity of (jA2l) , ODEs (lAll) and (IA3I) are obviously topologically 
equivalent. Moreover, they are also geometrically equivalent due to the discussion above. 


